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Abstract 

Wavefunctions of a heavy-light quark (Q, q) system described by 
a covariant Dirac hamiltonian are analyzed. By assuming that the 
confinement potential is a Lorentz scalar {S), the slope of the Isgur- 
Wise function is calculated at zero recoil point. The result obtained 
is £'(1) = —0.93 ± 0.05. This means that the solutions are perfectly 
consistent. If relativistic corrections in the light quark wave functions 
are included the result is £'(1) = —1.01 ± 0.04. From heavy-light 
data this suggests that if relativistic effects are considered, scalar con- 
finement is reliable in low orbital angular momentum states. 
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A covariant Dirac-like hamiltonian describing a heavy quark-light quark 
(Q, q) system has been introduced in Ref. [|]]. This hamiltonian contains all 
of the relativistic corrections. These effects come mainly from the recoil of the 
heavy quark. The model has been successful in applications involving both 
hydrogenlike (B-, D-) meson spectroscopy and also in calculations of weak 
mixing angles [[| from exclusive semileptonic B-decays. The hamiltonian 
of Ref. j l] contains two potentials, one of them is a central linear Lorentz 
scalar potential (S) which describes dynamically the confinement, and the 
other is a Coulomb potential (U) describing the short range color interquark 
interaction. To be strong enough, the non-perturbative scalar potential S 
does not feel any relativistic corrections but the Coulomb potential U still 
is susceptible to them. In fact, in the model of Ref. |TJ the potential U 
incorporates relativistic corrections to first order in tu/Mq, where m is the 
light quark mass and Mq the heavy quark mass. 

One question about the above description still remains to be studied and 
it is the concerning to the consistency of its solutions. For this reason, the 
purpose of this work is to examine the solutions of [|Xj] when the system (Q, q) 
is in low orbital angular momentum states. In order to do this, we subject the 
wavefunctions to several phenomenological tests. For instance, we check that 
they are normalizable. The way in which this is done is through an approach 
introduced in Ref. ||. Another thing which is checked in the present work 
is that the solutions account for the so called Klein paradox. This effect 
consists in the mixing of positive-negative energy states [|J. The third and 
last test employed on the wavefunctions is to calculate the slope of the Isgur - 
Wise (IW) function at zero recoil point, £'(1) 0. This quantity is calculated 
by assuming a scalar confinement potential (S) and also in the heavy quark 
symmetry limit (Mq — > oo). The value obtained for the slope is compared 
to other values of this quantity previously found using different approaches. 
After we check that the wavefunctions of Ref. ij are consistent, then we 
include relativistic corrections. These corrected solutions are employed to 
make the corrections to £'(1) that result from a finite Mq. Finally, the 
present work is concluded by discussing the nature of confinement in a (Q, q) 
system in low angular momentum states. 

Let us first check the normalizability of the solutions of Ref. [|J. To 
do this, we follow an approach similar to that of Ref. 0. In |J it was 
pointed out that a relativistic quarkonia system (q, q) with a vector confining 
potential V v = f3\f32K v r ( k v > ) could have unphysical non-normalizable 
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solutions. On the other hand, if the confining potential is of a Lorentz scalar 
nature, V s = K s r ( k s > ), this problem would not arise. If one uses both 
kind of potentials, it is necessary that scalar confinement be stronger than 
vector confinement (e.g. k s > k v ) in order to avoid singularities in the wave 
function. In a possible singularity in the (q, q) wave function is avoided 
through the use of a Salpeter (no-pair) equation. 

In order to illustrate the procedure, let us consider a simple scheme where 
the (Q, q) bound system is described by a Dirac equation. We assume that 
the only interaction existing between the quarks is that of confinement and 
that the confinement potential is composed of two parts. One of them is 
scalar, S = k s r (k s > 0), and the other vector, V = k v r (k s > 0). In this 
way, the respective equation in the CM. system is 



a 



p + m(3 + f3S + V ib = Etp, 



(1) 



where p is the light quark momentum. 
Let us decompose the wave function as 



ib = ip+ + = p±^, (2) 

where P ± = (1 ± (3)/ 2 are the standard projection operators into the upper 
and lower components of the wave function. Then, by using (1) and (2), we 
obtain 



ib~ = ^ — (3) 

V E + m + (k s - k v ) r v K } 

From the last expression it is evident that if E > and k v > k s , the lower 
sector of the wave function has a singularity at r = (E + m)/(k v — k s ). 
Consequently, the norm of the wave function 



< ib | ib >=< ib + | ib + > + < tb | tb >, (4) 

is not finite. In particular, we can observe from Eq. (3) that with a vector 
confinement (S = 0), the wavefunctions would not be finite. But, if the 
confinement is scalar (V = 0), the solutions are finite. 
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Let us turn now to examine the normalization of the wavefunctions of 
the covariant hamiltonian which was introduced in Ref. [ l] for describing a 
(Q, q) system. This hamiltonian is 



a ■ p + m(3 + 



P 



2M, 



+ M Q +U(r) + 



U(r) 



2M, 



Q 



a • p + a • f f • pj + f3S(r) 



iff = Eip, 



(5) 



where p is the light quark momentum, V s (r) = (3S(r) = f3n s r [k s > 0) is 
the scalar confining linear potential and U(r) = —^/r (£ > 0) is a color 



Coulomb-like potential. Note that terms proportional to 
this hamiltonian arise from the recoil of the heavy quark. 
Using Eqs. (2) and (5) we obtain 



U(r) 
2 M 



and to Tr^TT in 



2 Mr 



a ■ p + 2^ f a-p + a-rr-p j 

j — _ L 

E + m -2M^~ M Q- U(r) + S(r) 

where \E\ > Mq + m. Since —U(r) + V(r) > then for e = E — Mq — 
m > the wave function ip~ could have an unphysical singularity, which 
would depend strictly on the size of the term ^f^- Q- According to Heavy 
Quark Effective Theory (HQET) || this term represents the exchanged mo- 
mentum between the bound quarks. It acquires a maximum value in the 
nonperturbative regime (quarks and gluons inside the meson) ||. Thus, if 
we denote by R the radius of the hadron, then from HQET it follows that 



1 (6) 



1/R fms' 1 = 0.2/ R GeV. (7) 



2 Mq 

Thus, if we define the denominator of Eq. (6) as 



1 In the infinitly heavy quark mass limit this term does not exist. Consequently if e > 0, 
the solutions of Eq. (5) are perfectly well defined for any value of r. 
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h{r) = e + 2m - U(r) + S(r) - (8) 

and use (7) together with U = — £/r fms^ 1 = — 0.2£/ | r | GeV, 
S = K s r fins 1 = 5 k s \ r | Ge\^ _1 , and a; = r/R , we obtain 

/i(r) = e + 2m + (k s R 2 x 2 - x + £)/r. (9) 

In order to show that h(r) is positive for any value of r, let us note 
that the quantity between parenthesis in the last equation is positive at 
x = and it does not take negative values for any value of r providing 
that AR 2 k s £ > 1. For « = 0.2GW- 2 and f = 0.445 [§, this last 
condition means R > 1.6 GeV^ 1 = 0.32 fms. Since this condition is very 
reasonable for a hadronic radius, we can conclude that h(r) > for any 
value of r. In particular, we have shown that the denominator h(r) of Eq. 
(6) is not zero for any value of r. Therefore, the solutions of Eq. (5) are 
finite and normalizable. 

Let us proceed now with the second test of the solutions of the hamilto- 
nian introduced in Ref. [jXJ] . For this purpose we impose the s-wave classical 
turning points condition p = in Eq. (5). Once we do this, we obtain 
/3 ( m + S ) + Mq + U = E which becomes 

E + - M Q = m + S + U, (10) 



E_ - M Q = -m - S + U, (11) 

for the positive {(3 = 1) and negative ((3 = —1) energy states respectively^]. 

At this stage there are two possibilities, either U = or U ^ 0. If we 
turn off the Coulomb interaction ( U = ) in the above equations we are 

2 The condition p = can be seen from the MIT bag model point of view as the 
boundary condition that prevents light quark q current flux leaving through the meson 
bag. In this way, we expect typical values for the returning point (r r p ) of order r r p ' > 
R ~ 0.78. fm 0-@ 
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reproducing the same analysis as Ref. M. In this work the confinement in 
a Dirac equation was analyzed by assuming that the confinement potentials 
were of both kinds, scalar (S) and vector (V). It was found that the structure 
of the confinement must be scalar (V = 0) in order to avoid the Klein paradox 
(mixing of positive with negative energy states). For the above reason, in the 
present work we are just considering the second and more realistic situation 
where (7^0. Thus, if we turn on the color Coulomb interaction and plot 
Eqs. (10) and (11) with S(r) = n s r and U{r) = — £ / r for two different 
values of the light quark masses m = , and m = 0.5 GeV, we obtain 
the plot of Fig. 1. As can be seen from this figure, the singularity in the 
Coulomb potential at r = distorts and mixes the positive and negative 
energy states inside the domain of perturbative physics i. e. r <C R Q 
But (as should be expected of a good and consistent confinement potential) 
this effect dissappears in the confinement region r ~ R. In other words, 
Fig. 1 indicates that the scalar confining potential S of Eq. (5) accounts for 
the Klein paradox in the confining region r ~ R independent of the value 
of the light quark mass. Furthermore, we note also from Fig. 1 that the 
critical values of the radial coordinate (r c ) for which the physical condition 
E + — Mq — m > begins to be valid, independently of the value of m, 
are those such that r > r c where r c ~ 0.3 fms. If we think of r c as the 
critical value dividing the perturbative and the non perturbative region, it 
follows that the critical energy necessary to reach the non perturbative region 
(r > r c ) would be e c ~ l/r c ~ 0.67 GeV . Another sign of consistency is 
that precisely M D{Ds) — M c ~ M B — M b ~ 0.66 GeV ~ e c [§. 

It is easy to see that a confining linear vector potential (V = n v r) 
without a scalar potential (5* = 0) in Eq. (5) is automatically discarded. In 
fact, in this case the s-wave returning point condition p = would yield 

E + - M Q = m + V + U, (12) 

E_ - M Q = -m + V + U. (13) 

For large values of r say r ~ R where the Coulomb potential vanishes 
{U ~ 0), then Eqs. (12) and (13) become 

3 In fact, as is well known, this mixing inside the confinement region r <!C R is respon- 
sible for the bound state spectra in a hydrogenlike system (Q,q). 
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E_ ~ E + ~ M Q + V > ; 



(14) 



This result is unphysical since it means that the negative and positive energy 
states are mixed in the non-perturbative region. In this way, we can conclude 
that only the scalar potential of Eq. (5) can account for the Klein paradox. 
This result makes the solutions of (5) consistent. 

In order to continue the analysis of the solutions of Eq. (5) let us turn 
now to the calculation of the slope of the IW function at zero recoil point. 

By solving Eq. (5) in the heavy quark symmetry limit (Mq — ► oo), we 
find the light quark wavef unctions. With these solutions the IW function is 
calculated through the expression [|] 

fa) = ZJTlOo(2£ g V^)>> ( 15 ) 
where E q is the light quark energy. The average is defined as 

(/> = f n drr^(ryf(r)^(r), (16) 

where Q is the spatial region explored by the wave function. In this work 
Q is taken as a sphere of radius R . Likewise, for the calculation of the 
slope of the IW function at zero recoil point, we use both the solutions of (5) 
in the heavy quark symmetry limit and the expression || 

£'(1) = -(§ + |^ 2 <r 2 >). (17) 
Once we do the above we obtain 

f(l) ~ -0.93 ± 0.05. (18) 

We have allowed the light quark mass to run from to 0.25 GeV and the 
hadron radius as R = 0.78 fm = 3.9 GeV- 1 [Q. 

The above value for £,'(!) is consistent with values of this quantity previ- 
ously obtained in other works. For instance, in Ref. [[| £'(1) — — 0-90 was 
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obtained by solving the Dirac equation with scalar confinement (S) through 
a variational method. While in Ref. || by analyzing sum rules it was found 
£'(1) — — 0.65. In [|HJ a relativistic flux tube mesonic model was employed 
and gave the result £'(1) — — 0.93, which is exactly the same value found 
in the present work. The coincidences between (18) with values for £'(1) 
obtained with such different approaches allow us to assert that solutions of 
(5) can go through this last test. 

Since the solutions of (5) have passed satisfactorily the three different 
tests above cited we can conclude that the model introduced in Ref. U is 
perfectly consistent to describe a (Q,q) system in low orbital angular mo- 
mentum states. 

Now that we have checked the consistency of Eq. (5), let us employ its 
solutions to discuss the nature of confinement in a (Q, q) system. 

As was pointed out in Ref. ||, the value of the IW function constitutes a 
sensitive test for the kind of confinement potential. According to this work, 
there is an apparent conflict at the moment between the values for £'(1) 
obtained from a Dirac equation with confining scalar on the one hand and 
that extracted from heavy-light 'data' on the other. The values of £'(1) from 
the data are higher than that calculated with a Dirac equation. To quote 
just a few values obtained from different Lattice QCD (LQCD) calculations, 
S'Mlqcd = -1.0 Q, ?{1)lqcd = -1.2 0, and ^'{1) LQ cd = -1.16. 
Such a large discrepancy between £'(1) and £'(1)lqcd in || questions the 
scheme of a scalar confining potential in a Dirac equation to describe a (Q, q) 
system. In fact, the authors of 0] propose that if instead of a Dirac equation 
with scalar confining potential a Salpeter (no-pair) equation with a vector 
confining potential is used to calculate £'(1), the values obtained for £'(1) 
would be in better agreement with heavy-light 'data'. The respective values 
for £'(1) found in j| with these two different models are ^'{1) Dirac — —-90 
and no-pair — —1-2. From these values they then argue in favor of a 
no-pair equation with vector confinement. 

We want to stress at this point that the values obtained in |4j] for £'(1) 
do not contain in them corrections coming from assuming Mq finite. Fur- 
thermore, it is claimed in this work that if one includes these corrections in 
the light quark wave functions, the scheme of a Dirac equation with vector 
confining potential would be the correct one for the description of a (Q, q) 
system. That these corrections are necessary is evident from the following 
three facts: 1. The limit Mq — > oo where both the IW function and its slope 
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are properly defined are just a good approximation to reality. 2. Even in 
the limit Mq — > oo, the IW function and its slope are of a non-perturbative 
nature. Consequently their values rely on parameters (quark masses, Bag 
radius, tensions, etc) independently of the approach employed. 3. A good 
model for the description of a (Q, q) is that which retains its validity in the 
limit Mq — ► oo. Furthermore, this model can be seen as a good heavy quark 
picture modified by finite relativistic corrections. 

As we stated above, we consider it important to include relativistic cor- 
rections to the slope of the IW form factor at zero recoil point. The way 
we incorporate these relativistic corrections to £'(1) is as follows. First, we 
obtain the light quark wavefunction from Eq. (5) with finite Mq using for 
this purpose the same parameters as in Ref. [|IJ. These solutions are then 
substituted in Eq. (17). Once the above procedure is implemented, we obtain 

£,'{±)corr = -1.01 ± 0.04, (19) 

where the light quark mass has taken the same range of values as in Eq. (18). 

As can be seen in Eq. (19) the value for £'(1) with relativistic corrections 
is in considerably better agreement with heavy-light data than those without 
relativistic corrections in Eq. (18) and Ref. @]. The corrected value for 
£'(1) of Eq. (19) modifies sustantially any conclusion about the nature of 
confinement. In fact, this result induces us to conclude that a Dirac equation 
with scalar confinement constitutes a good way of describing a (Q,q) system 
in low orbital angular momentum states providing relativistic corrections are 
taken into account. 
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FIGURES 



Figure 1. Behavior of the right hand side of Eqs. (10) and (11) with U(r) = 
— £/r and S(r) = K s r for two different light quark masses: la), m = and 
lb), m = 0.5 GeV. The values employed for £ and the same of Ref. 

[§ e. g. i = 0.445 and k s = 0.2 GeV~ 2 . 
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